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对 流 一 扩散 一 反应 方程 的 变 分 多 尺度 解法 * 


朱 海 涛 ， 欧阳 洁 
(西北 工业 大 学 应 用 数学 系 ， 西 安 710072) 


搞 要 : 根据 变 分 多 尺度 的 思想 求解 了 对 流 项 和 反应 项 占 优 的 对 流 一 扩散 一 反应 方程 。 在 变 分 多 尺度 思想 
的 理论 框架 内 ， 推 导 了 附加 于 Galerkin 变 分 弱 形式 的 稳定 化 结构 和 具体 的 稳定 化 系数 ， 培 述 了 
这 种 稳定 化 结构 和 经 典 的 SUPG 稳定 化 结构 之 间 的 关系 ; 数值 算 例 表明 ， 该 稳定 化 系数 可 以 适 
应 均匀 和 非 均匀 的 计算 网 格 。 通 过 网 格 的 恰当 加 密 ， 变 分 多 尺度 方法 消除 了 算 例 中 的 数值 伪 振 
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1 引言 


用 标准 Galerkin 方 法 求解 对 流 一 扩散 一 反应 方程 ， 当 对 流 项 或 反应 项 占 优 时 ， 数 值 结 
果 会 伴随 着 剧烈 的 数值 伪 振 荡 。 通 过 在 Galerkin 变 分 弱 形式 的 左 端 加 上 一 项 基于 原 偏 微 
分 方程 残 差 的 稳定 项 而 产生 的 稳定 化 类 方法 ， 既 使 该 类 数值 计算 方法 获得 了 相 容 性 ， 又 
获得 了 稳定 性 四 ， 从 而 保证 了 稳定 化 类 方法 的 收敛 性 。 这 类 方法 包括 : SUPG (Streamline 
Upwind /Patrov-Galerkin)、GLS (Galerkin Least Squares). USFEM (Unusual Stabilized Finite 
Element Methods) 等 。 但 是 ， 这 类 方法 的 缺点 是 最 优 的 稳定 化 系数 很 难 确定 。 上 世纪 90 年 
{X, Hughes 考虑 了 各 种 稳定 化 方法 ， 提 出 了 变 分 多 尺度 方法 的 基本 理论 框架 3 站。 本 文 根 据 变 
分 多 尺度 的 思想 求解 了 对 流 项 和 反应 项 占 优 的 对 流 一 扩散 一 反应 方程 ， 主 要 的 工作 是 : 1) 在 
变 分 多 尺度 的 理论 框架 内 ， 用 泡 函 数 近似 “ 细 ” 尺 度 上 的 解 ， 推 导 了 用 变 分 多 尺度 方法 求解 该 
方程 的 附加 于 Galerkin 变 分 弱 形 式 的 稳定 化 结构 及 对 应 的 稳定 化 系数 ，2) 阐述 了 在 线性 和 双 
线性 的 情形 下 ， 变 分 多 尺度 方法 的 稳定 化 结构 和 经 典 的 SUPG 稳定 化 结构 之 间 的 关系 ; 3) HE 
于 变 分 多 尺度 思想 求解 了 两 个 具有 代表 性 的 数值 算 例 ， 并 通过 网 格 加 密 技术 与 变 分 多 尺度 方法 
的 耦合 ， 消 除了 算 例 中 由 对 流 项 和 反应 项 占 优 引 起 的 数值 伪 振荡 。 


2 ”对 流 -扩散 -反应 方程 
考虑 如 下 的 对 流 一 扩散 一 反应 方程 


a-Vu—kAut+su=f EQ 内 ， 
u=0 在 6Q E, 
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其 中 a 为 已 知 的 对 流速 度 ; 上 是 扩散 系数 ，k > 0; s 是 已 知 的 反应 项 系数 ， 了 是 已 知 的 源 项 。 
记 为 求解 区 域 Q 的 下 规章 分 ，K € G4 是 C4 中 的 任 -有 限 单元 ，hx Roe K KH 


h= max {hi}, Vi = {un € HG(Q) | vale € P'(K), VK € Ch}, 


PK) 是 单元 K 上 连续 的 线性 或 双 线 性 多 项 式 函数 空间 。 
方程 (1) 对 应 的 标准 Galerkin 离散 格式 为 : Ru, € Vi C HEN), WE 


(a: Vup,vVh)o 十 k(Vun, Vuro 十 (sun, vala = (fun), Vur EVIC Hi, (2) 


这 里 (u,v) = /5 uvdz， 除 非特 别 声 明 ，D = Q。 当 kk < |a Rk < s 时 ， 上 述 求 解 格式 会 伴 
随 着 剧烈 的 无 任何 物理 意义 的 数值 伪 振 荡 。 本 文 根 据 变 分 多 尺度 的 基本 思想 ， 用 泡 函 数 近 似 
“ 细 ” 尺 度 上 的 解 ， 引 入 了 类 似 经 典 SUPG 稳定 化 方法 的 稳定 项 ， 增 加 数值 求解 格式 的 稳定 
性 。 


3 ”对 流 -扩散 -反应 方程 的 变 分 多 尺度 分 析 


iT =Uke6OK, B= {b| b € H}/Vi, b=0, ÆT E}, M = Ukec (intK)。 在 变 分 多 尺 

度 的 理论 框架 内 ， 首 先 把 试验 函数 和 检验 函数 分 解 为 “ 粗 、 细 ”两 种 尺度 ， 即 
u(t) = a(z) + wa), (3) 
v(x) = U(x) + v(x), (4) 


HP al) Mo(c) 分 别 为 “ 粗 ” 尺 度 上 的 试验 函数 和 检验 函数 ，w(z) 和 w (zx) 分 别 为 “ 细 ” 尺 
度 上 的 试验 函数 和 检验 函数 。ia(z) E Vi, 5(z) € Vi; w(x) € B, v(x) € B, At, u(x) € 
Vi @ B， 即 我 们 用 BB 中 的 函数 近似 “ 细 ” 尺 度 上 的 解 ， 后 面 将 会 看 到 ， 所 谓 的 “ 细 ” 尺 度 的 解 
其 实 是 “ 粗 ” 尺 度 解 的 误差 。 把 式 (3)、(4) 代入 式 (2) 中 可 得 


(a-V(itw),o+v') +k(V(@ +u), VT+v)) + (s(a +u), +v) =(f,0+0'), (5) 

HV, Al B WREEKT, HAER (5) 分 解 为 “ 粗 、 细 ”两 种 尺度 ， 即 
(a: V(a +u), 8) +k(V(a+u’), Vo) + (s(ä +u’), a) = (f,0), (6) 
(a: V(@ +u), v) +k(VG+w’), Vv’) + (s(@ +u’), v’) = (f,0’). (7) 


31 细 尺 度 上 方程 的 求解 
Fu’ OK 上 为 零 ， 式 (7) 可 以 写 为 如 下 的 形式 


(a-Vitu'),v')o, +k(VG+u'), Vo) o + (s(@+v),v')o, = (w, (8) 
其 中 
(u,v) = uvdz. 
woe Ie 
A (8) 等 价 于 


(a - Vi, vo + (a: Vu’, v) + k(Vū, Vu a 


+ k(Vu', Vv’) + (st, var + (su, v) = (fvw, (9) 
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对 上 式 左 端 第 三 项 应 用 散 度 定理 ， 整 理 可 得 
(a-Vu',v' jw +k(Vu', Vv Ja + (su, vq = —(a- Va kAL sū — f, vw, (10) 
R (10) 体现 了 变 分 多 尺度 方法 的 基本 思想 ， 该 式 等 价 于 
(a: Vu — kAv + su’)q = —(a- Va — kAū + st — flo, (11) 


即 ，“ 细 ”尺度 解 是 “ 粗 ” 尺 度 解 的 误差 ，“ 粗 、 细 ”尺度 上 的 解 的 和 为 原 偏 微分 方程 的 精确 
解 。 需 要 指出 的 是 : 式 (11) 只 在 单元 内 部 成 立 。 
wWe=a-Vu-kAt+su-—f, aR “fA” RES RN RAT ENR. (10) 即 为 


(a: Vu’, vw + k(Vu', Vv’) + (su’, v) = 一 (su or， (12) 


上 式 就 是 我 们 要 求解 的 细 尺 度 上 的 变 分 问题 ， 它 与 式 (11) 等 价 。 求 解 式 (12) 有 多 种 方法 ， 
如 Green 函数 法 2,3， 双 层 有 限 元 方法 四 等 ， 本 文采 用 泡 函 数 近似 “ 细 ” 尺 度 解 求解 式 (12)。 
与 稳定 化 方法 相对 应 ， 变 分 多 尺度 方法 在 “ 细 ” 尺 度 上 的 解 将 最 终 影 响 稳 定 化 因子 的 确定 。 为 
了 在 稳定 化 因子 中 考虑 对 流 项 的 影响 ， 求 解 式 (12) 的 过 程 中 选取 了 不 同 的 试验 函数 和 检验 函 
数 。 
为 简单 起 见 ， 假 设 试验 函数 和 检验 函数 均 只 用 一 个 基 函 数 ， 分 别 为 站 和 上 既 ， 则 试验 函数 和 
检验 函数 分 别 为 : ul = ubi, v = vpb§, RAR (12) 可 得 


up[(a- Vbi, b5)ar + k(Vbi, VoS)a + (sbi, b3)ar] ve = —(e, bs)ov, (13) 
由 w 的 任意 性 ， 有 
up = —(€, b$) [(a : VOF, bJa + b(Vb§, Vo§)or + (sbi, bar], (14) 


即 
u = 好 { = (e, b5) [(a Voi, b5) 二 k(Vbĵ, Voz) + (5 好 ， bi) ro (15) 


为 了 引入 基于 原 偏 微分 方程 残 差 的 稳定 项 ， 保 证 整体 求解 格式 的 相 容 性 ， 用 常数 近似 由 “ 粗 ” 
尺度 上 的 解 引 起 的 原 偏 微分 方程 在 每 个 单元 内 的 残 差 ， 引 入 记号 r，“ 细 ”尺度 上 的 解 可 简化 
为 


ul = —T(a: Vi — kAu + st — f), (16) 
T = 0{(1,08)% [(a- Vo, b5) + k(VB$, Vga + (sbi, §)a"] “}. (17) 


3.2 ” 粗 尺 度 上 方程 的 求解 
考虑 式 (6)“ 粗 ”尺度 方程 的 求解 。 式 (6) 等 价 于 


(a- Va, v) + (a- Vu', 可 十 有 ViVi 十 ECVwV 可 十 (5 一 可 十 (su = (f, 0), (18) 
对 上 式 左 端 第 二 、 四 项 应 用 散 度 定理 ， 整 理 可 得 


(a- Va, T) + k(Vd, Vv) + (st, 6) + (u’, -a- Vi — KAD + sv) = (f,50), (19) 
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把 式 (16),(17) 代入 式 (19) 中 ， 有 
(a- Vi, U) + k(Vu, Vo) + (st, 0) 
+(-r(a- Vū-— kAūŭ + su — f),-a- Vi— kAd+ sd) = (f,0). (20) 


3.3 稳定 化 结构 和 稳定 化 因子 
稳定 化 类 方法 的 计算 格式 由 标准 的 Galerkin 离散 格式 和 稳定 项 两 部 分 组 成 ;稳定 项 由 残 差 
项 、 稳 定 化 因子 和 稳定 化 函数 组 成 ， 即 一 般 的 稳定 化 计算 格式 为 


G(u,v) + >》 (el(u), 7P(v)) = (fv). (21) 
K 


其 中 G(u, v) 为 标准 的 Galerkin RAH, elu) 为 原 偏 微分 方程 的 残 差 ，P(w) 为 稳定 化 函数 ，7 为 
稳定 化 系数 。 考 虑 经 典 的 稳定 化 方法 SUPG， 有 


P(v) =a- Vv, (22) 


7 是 一 个 根据 经 验 人 为 设计 的 一 个 稳定 化 因子 ， 有 具体 计算 公式 参考 文 [1,5]。SUPG 方法 引入 了 
基于 原 仿 微 分 方程 残 差 的 稳定 项 ， 因 此 可 以 保证 稳定 化 后 计算 格式 的 相 容 性 ， 关 于 SUPG 稳定 
化 方法 稳定 性 的 证 明 可 参见 文 15]， 由 相 容 性 和 稳定 性 保证 了 数值 求解 格式 的 收敛 性 。 

考虑 变 分 多 尺度 方法 引入 的 稳定 项 ， 由 式 (20) 知 


P(v) =a-v+kAv — sv, (23) 
7 由 式 (17) 给 定 。 在 实际 计算 中 ， 出 于 简单 性 的 考虑 我 们 忽略 了 第 三 项 的 影响 ， 即 取 
P(v) = a:v+ kAyv, (24) 


显然 ， 在 线性 和 双 线 性 的 情形 下 ， 变 分 多 尺度 方法 引入 的 稳定 项 和 SUPG 稳定 化 方法 的 稳定 项 
有 相同 的 结构 。 由 于 该 方法 引入 了 基于 原 偏 微分 方程 残 差 的 稳定 项 ， 因 此 ， 可 以 保证 求解 格式 
的 相 容 性 。 

与 传统 稳定 化 方法 的 稳定 化 因子 不 同 ， 变 分 多 尺度 方法 引入 的 稳定 化 因子 不 是 一 个 人 为 设 
计 的 常数 ， 而 是 一 个 与 “ 细 ” 尺 度 解 密切 相关 的 函数 。 当 对 流 项 占 优 时 ，7t 是 O(hjlal); 当 扩 
散 项 占 优 时 ，7 是 O(h2/k); 当 反 应 项 占 优 时 ，r 是 O(1/s)， 这 与 文献 [47] 中 的 分 析 一 致 。 

基于 上 述 分 析 ， 采 用 线性 和 双 线 性 有 限 单 元 求解 方程 (1) 的 变 分 多 尺度 公式 为 


(a: Va, o) + k(VG, Vo) + (sù, 5) 


+ (r(a- Va — kAū + sū),a - Vū + kAd) = (f, 0) + (f, Tr(a- Vū + kAD)). (25) 


4 ”数值 算 例 


对 流 一 扩散 一 反应 方程 不 仅 在 边界 上 可 以 形成 尺度 很 小 的 边界 层 效 应 ， 而 且 当 边界 条 件 存 
在 间断 时 还 会 出 现 尺度 非常 小 的 内 部 层 效 应 。 本 部 分 我 们 考虑 两 个 具有 代表 性 的 数值 算 例 ， 第 
一 个 算 例 为 连续 的 边界 条 件 ， 主 要 形成 边界 层 效 应 ;第 二 个 算 例 为 间断 的 边界 条 件 ， 形 成 边界 
层 效应 和 内 部 层 效 应 。 

4.1 数值 算 例 1 
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首先 考虑 如 下 的 数值 算 例 


a-Vu-—kAu+su=1 f Q= [0,1] x [0,1] A, 
u=0 在 09 上， 


Hk = 1074, a= |al(cos(z/3),sin(7/3)), |al R s 分 别 取 以 下 三 种 情形 : 

(a) |al = 1.0, s = 0.0001， 对 应 于 对 流 占 优 时 的 情形 ; 

(b) ja| = 0.5，s = 1.0， 对 应 于 对 流 、 反 应 同时 占 优 ; 

(c) Ja] = 0.0001,，s = 1.0， 对 应 于 反应 占 优 时 的 情形 ; 

图 1， 图 2 分 别 是 经 典 SUPG 稳定 化 方法 和 变 分 多 尺度 方法 的 计算 结果 ， 其 中 计算 网 格 
为 20 x 20 的 均匀 矩形 网 格 。 在 (alj)、(b)、(c) 三 种 情形 下 ， 标 准 的 Galerkin 方 法 均 存 在 严重 
的 数值 念 振荡。 图 1 是 对 应 于 (a) (b) (c) 三 种 情形 的 经 典 SUPG 稳定 化 方法 的 数值 结果 ; 
图 2 是 对 应 于 (a)、(b)、(c) 三 种 情形 的 变 分 多 尺度 (Variational Multiscale, 简称 VM) 方法 的 数 
值 结 果 。 可 以 看 出 ， 当 对 流 占 优 时 ， 在 边界 处 变 分 多 尺度 方法 有 较 SUPG 更 剧烈 的 数值 伪 振 
荡 ， 原 因 是 实际 计算 中 采用 的 简化 后 的 “ 细 ” 尺 度 模 型 在 较 粗 的 计算 网 格 上 没有 很 好 地 捕捉 到 
具有 极 细 尺度 特征 的 边界 层 效 应 的 物理 信息 。 
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(a) lal = 1.0, s = 0.0001 (b) ja| =0.5, s=1.0 (c) lal = 0.0001, s=1.0 


图 1: 均匀 网 格 上 SUPG 方法 的 数值 解 
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(a) lal = 1.0, s = 0.0001 (b) lal=0.5, s= 1.0 (c) lal = 0.0001, s = 1.0 


图 2: 均匀 网 格 上 VM 方法 的 数值 解 
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为 了 消除 边界 上 存在 的 数值 伪 振 荡 ， 我 们 设计 了 在 边界 处 加 密 的 非 均 匀 网 格 ， 采 用 的 20 x 
20 的 非 均匀 和 矩形 计算 网 格 如 图 3 所 示 。 图 4 是 对 应 于 (a)、(b)、(c) 三 种 情形 的 经 典 SUPG 稳定 
化 方法 的 数值 结果 ; 图 5 是 对 应 于 (a)、(b)、(c) 三 种 情形 的 变 分 多 尺度 方法 的 数值 结果 。 由 
图 5 可 见 ， 我 们 得 到 的 结果 明显 优 于 文献 [6] 的 结果 ， 后 者 使 用 SUPG 稳定 化 方法 ， 尽 管 采用 
了 42 x 42 的 非 均匀 计算 网 格 ， 但 仍然 没有 消除 在 点 (1.0, 1.0) 处 的 数值 伪 振 荡 ， 而 图 5 的 结果 
完全 地 消除 了 边界 处 的 数值 伪 振 荡 。 与 此 同时 ， 对 比 图 4 和 图 5 可 发 现 ， 与 变 分 多 尺度 解 相 
比 ，SUPG 解 的 精度 较 低 。 其 原因 是 SUPG 解 的 精度 与 稳定 化 因子 密切 相关 ， 而 稳定 化 因子 又 
与 网 格 尺寸 、 网 格 疏 密 有 紧密 的 联系 。 人 为 设计 的 稳定 化 因子 对 非 均匀 的 计算 网 格 缺乏 弹性 ; 
而 由 变 分 多 尺度 方法 计算 出 来 的 稳定 化 因子 实际 是 一 个 与 “ 细 ” 尺 度 上 解密 切 相关 的 函数 ， 因 
此 可 以 很 好 地 适应 均匀 和 非 均匀 网 格 。 “ 


图 3: 算 例 1 的 非 均匀 计算 网 格 


(b) la| =0.5, s =1.0 (c) Jal = 0.0001, s = 1.0 


图 4: 非 均匀 网 格 上 SUPG 方法 的 数值 解 


(a) lal = 10, s = 0.0001 (b) lal = 0.5, s= 1.0 (c) Jal = 0.0001, s = 1.0 


图 5: 非 均匀 网 格 上 VM 方法 的 数值 解 
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4.2 数值 算 例 2 
考虑 如 下 具有 间断 边界 条 件 的 数值 算 例 


a-Vu-—kAu+su=0 # Q= [0,1] x [0,1] A, 
u=1 Æ c=0, X y=1 E, 
u=0 在 z=1， X y=0 E, 


HH k = 10-4, a= |al(cos(r/4),sin(r/4))，|al 和 s 分 别 取 以 下 三 种 情形 ， 

(a) la| = 1.0，s = 0.0001， 对 应 于 对 流 占 优 时 的 情形 ; 

(b) |a| = 0.5，s = 1.0， 对 应 于 对 流 、 反 应 同时 占 优 ; 

(c) |a| = 0.0001, s = 1.0， 对 应 于 反应 占 优 时 的 情形 ; 

在 (a)，(b)，(c) 三 种 情形 下 ， 标 准 的 Galerkin 方法 均 伴 随 着 剧烈 的 数值 伪 振 荡 ， 图 6 是 对 
应 于 以 上 (a)，(b)，(c) 三 种 情形 的 变 分 多 尺度 方法 的 计算 结果 ， 其 中 情形 (a) 主要 形成 内 部 层 
效应 ， 因 此 采用 40 x 40 的 均匀 网 格 ， 情 形 (b) 和 (c) 均 采 用 边界 附近 局 部 加 密 的 40 x 40 非 均匀 
网 格 。 可 以 看 出 ， 变 分 多 尺度 方法 消除 了 由 边界 层 效应 和 内 部 层 效 应 引起 的 数值 伪 振 荡 。 


(a) |a| = 1.0, s = 0.0001 (b) la] = 0.5, s = 1.0 (c) |a| = 0.0001, s = 1.0 


图 6: 算 例 2VM 方法 的 数值 解 


5 ”结论 


变 分 多 尺度 方法 是 一 种 有 效 处 理 具有 多 尺度 效应 的 数学 物理 问题 的 数值 方法 。 在 变 分 多 尺 
度 方法 的 框架 内 可 以 推导 出 经 典 的 稳定 化 方法 ， 如 SUPG、GLS 稳定 化 方法 等 ， 与 此 同时 变 分 
多 尺度 方法 引入 了 和 经 典 稳定 化 方法 结构 相同 的 稳定 项 ， 保 证 了 计算 格式 与 原 偏 微分 方程 的 相 
容 性 。 除 此 之 外 ， 在 该 框架 内 可 以 自动 地 计算 出 稳定 化 因子 ， 该 稳定 化 因子 实际 上 是 一 个 与 
“ 细 ” 尺 度 上 的 解密 切 相关 的 函数 ， 它 较 人 为 根据 经 验 设计 的 稳定 化 系数 具有 更 强 的 适应 性 ， 
不 仅 适 用 于 均匀 网 格 ， 也 适应 非 均匀 网 格 。 通 过 网 格 加 密 技 术 与 变 分 多 尺度 方法 的 耦合 ， 可 以 
消除 在 对 流 -扩散 -反应 方程 中 由 边界 层 效 应 和 内 部 层 效 应 引起 的 数值 伪 振 荡 。 
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Variational Multiscale Method for the Advection-diffusion-reaction 


Equation 
ZHU Hai-tao, OUYANG Jie 


(Department of Applied Mathematics, Northwestern Polytechnical University, Xi’an 710072) 


Abstract: The advection-diffusion-reaction equation dominated by either advection or reaction is 


solved by using the variational multiscale method, which preserves the superior properties of the clas- 


sical stabilized methods such as SUPG, GLS. The variational multiscale method decomposes the scalar 


field into coarse scale and fine scale, respectively; then the fine scale incorporates into the coarse scale 


through a bubble function to correct the lack of stability of the coarse scale; at last, the coarse scale is 


solved by the finite element method. A distinctive feature of the proposed method is that the stabilized 


parameter can be calculated automatically which adapts to uniform and nonuniform grids. The spuri- 


ous oscillation in our numerical examples is eliminated by coupling the variational multiscale method 


with nonuniform grids. 
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